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abstract 



We predict branching ratios and CP asymmetries of the B — > KK decays using perturbative QCD factorization 
theorem, in which tree, penguin, and annihilation contributions, including both factorizable and nonfactorizable ones, 
are expressed as convolutions of hard six-quark amplitudes with universal meson wave functions. The unitarity angle 
(f>3 — 90° and the B and K meson wave functions extracted from experimental data of the B — > Ktt and irn decays are 
employed. Since the B — > KK decays are sensitive to final-state-interaction effects, the comparision of our predictions 
with future data can test the neglect of these effects in the above formalism. The CP asymmetry in the B^ — > K^K 
modes and the B% — > K ± K T branching ratios depend on annihilation and nonfactorizable amplitudes. The B — > KK 
*^ \ data can also verify the evaluation of these contributions. 

P-i 

I. INTRODUCTION 

The conventional approach to exclusive nonleptonic B meson decays relies on the factorization assumption (FA) , 
under which nonfactorizable and annihilation contributions are neglected and final-state-interaction (FSI) effects are 
assumed to be absent. Factorizable contributions are expressed as products of Wilson coefficients, meson decay con- 
stants, and hadronic transition form factors. Though analyses are simpler under this assumption, estimations of many 
important ingredients, such as tree and penguin (including electroweak penguin) contributions, and strong phases are 
not reliable. Moreover, the above naive FA suffers the problems of scale, infrared-cutoff and gauge dependences [||. 
It is also difficult to explain the observed branching ratios of the B — > J/ipK^ decays in the FA approach, to which 
nonfactorizable and factorizable contributions are of the same order || . 

Perturbative QCD (PQCD) factorization theorem for exclusive heavy-meson decays has been developed some time 
ago which goes beyond FA. PQCD is a method to separate hard components from a QCD process, which are 

treated by perturbation theory. Nonperturbative components are organized in the form of hadron wave functions, 
which can be extracted from experimental data. This prescription removes the infrared-cutoff dependence in PQCD. 
Since nonperturbative dynamics has been absorbed into wave functions, external quarks involved in hard amplitudes 
are on-shell, and gauge invariance of PQCD predictions is guaranteed. Contributions to hard parts from various 
topologies, such as tree, penguin and annihilation, including both factorizable and nonfactorizable contributions, can 
all be calculated. Without assuming FA, it is easy to achieve the scale independence in the PQCD approach. 

Despite of the above merrits of PQCD, an important subject, final state interaction (FSI), remains unsettled, 
which is nonperturbative but not universal. FSI effects in two-body decays have been assumed to be small. Though 
arguments and indications for this assumption have been supplied in fit], experimental justification is necessary. In 
this paper we shall propose to explore FSI effects by studying B — > KK decays. It will be explained in Sec. II that 
these decays are more sensitive to FSI effects compared to 5 Ktt and tttt decays. Employing the meson wave 
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functions and the unitarity angle 03 = 90° determined in J7j , we predict the branching ratios and the CP asymmetries 
of the B^ — > K^K , B^ — > K^K T and B^ — > K°K° modes. The comparision of our predictions with future data 
can be used to estimate the importance of FSI effects. In particular, large observed B® — > K ± K T branching ratios 
and CP asymmetry in the B^ — > K°K° modes will imply strong FSI effects. 

An essential difference between the FA and PQCD approaches is that annihilation and nonfactorizable amplitudes 
are neglected in the former, but calculable in the latter. It has been shown that annihilation contributions from 
the operator O5.6 with the (V ~ A)(V + A) structure, bypassing helicity suppression, are not negligible 0. These 
contributions, being mainly imaginary, result in CP asymmetries in the B — > tttt decays, which are much larger than 
those predicted in FA (8|^]. Hence, measurements of CP asymmetries will distinguish the two approaches ||. The 
B^ — > K^K a modes contain both annihilation amplitudes from 0$$ and nonfactorizable annihilation amplitudes 
from O12, such that they exhibit substantial CP asymmetry. The branching ratios of the B^ — > K^K T modes, 
involving only nonfactorizable annihilation amplitudes, can not be estimated, or are vanishingly small in FA. The 
data of these two decays can verify PQCD evaluation of annihilation and nonfactorizable contributions. 

FSI effects in the B — > tttt, Ktt and KK decays are compared in Sec. II. The PQCD formalism for annihilation and 
nonfactorizable contributions is reviewed in Sec. III. We present the factorization formulas of all the B — » KK modes 
in Sec. IV, and perform a numerical analysis in Sec. V. Section VI is the conclusion. 



II. FINAL STATE INTERACTION 



FSI is a subtle and complicated subject. Most estimates of FSI effects in the literature jlOj suffer ambiguities or 
difficulties. Kamal has pointed out that the enhancement of CP asymmetry in the B^ — » K^tt^ modes from order 
0.5 % up to order (10-20)% |l|] is due to an overestimation of FSI effects by a factor of 20jO|. The smallness of FSI 
effects has been put forward by Bjorken fl3| l based on the color-transparency argument |14| . The renormalization- 
group (RG) analysis of soft gluon exchanges among initial- and final-state mesons [Gl| has also indicated that FSI 
effects are not important in two-body B meson decays. These discussions have led us to ignore FSI effects in the 
PQCD formalism. For example, the charge exchange in the rescattering B + — ► K + tt° — ► K°tt + , regarded as occuring 
through short-distance quark-pair annihilation, is of higher order [Q. 

As stated in the Introduction, the neglect of FSI effects requires experimental justification. For this purpose, we 
propose to investigate the B — > KK decays, which are more sensitive to FSI effects compared with the B — > Ktt and 
tttt decays. Similar proposals have been presented in the literature [T^JT^ ] within the framework of SU (3) symmetry. 
We make our argument explicit by means of the general expression for the B — > tttt, Ktt and KK decay amplitudes, 



A = V U T + V U P U + V C P C + V t Pt • 



(1) 



The factors V q — V q d,V* b , q — u, c, and t, are the products of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements, 
T denotes the tree amplitude, and P q denote the penguin amplitudes arising from internal q-quark contributions. FSI 
effects have been included in the amplitudes T and P u ,c,t- 

Using the unitarity relation V c = — V u — Vti Eq. (0) is rewritten as 

A = V U (T + P U - P c ) + V t (P t - P c ) , 

V t P t - Pc 



V u (T + P u - P c ) 
V u (T + P u - P c ) 



1 



Vu T+Pu-Pc 



V u 



r(KK)t 



(2) 



A = V t {Pt - Pc) 
= V (P t - P c ) 



V u T + P u - P c 



V t 

Vu 

v t 



Pt-Pc 



(3) 



where 1Z are the ratios of different amplitudes and 6 the CP-conserving strong phases. 

Without FSI, the various amplitudes T and P u , c ,t, namely, the ratios 1Z and the strong phases 5 are calculable in 
PQCD. If FSI effects are important, they may change branching ratios or induce CP asymmetries of two-body B 
meson decays by varying R and S. For the B — > tttt decays, the ratio Vt/V u is of order unity, but R^n is small because 
of the large Wilson coefficients 01 = 6*2 + C\/N c in T, iV c being number of colors. The exception is the B^ — ■> 7r°7r° 
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mode, whose tree amplitude is proportional to the small Wilson coefficient 02 = C\ + CijN c . The ratio Rk-k may 
be large, but its coefficient V u /Vt ~ i?bA 2 , Rb and A being the Wolfenstein parameters defined in Sec. IV, is small. 
Therefore, FSI effects in the B — > irir and Kir decays are suppressed by \ ja\ and V u /Vt, respectively. On the other 
hand, the B — > 7T7t and AT7r decays have branching ratios of order I0 -5 , which are larger than those of the B — > KK 
decays (of order 10 -6 as calculated in Sec. V). It has been also-predicted in PQCD that the CP asymmetries in the 
B — > irir and Kir decays are large: 30 — 40% in the former |8||£j and 10 — 15% in the latter |?]]. These large values 
render FSI effects relatively mild. 

For the B — > KK decays, T arises only from small nonfactorizable annihilation diagrams for the B ± -> K±K n 
and B d — > K^K T modes, and vanishes for the B d — > K°K° modes. Furthermore, Vt/V u is of order unity. Hence, 
there is no suppression from the Wilson coefficients and from the CKM matrix elements, and FSI effects will be more 
significant. In the PQCD approach TZkk is close to unity, corresponding to the branching ratios of order 10 -6 for 
B ± K±K° and B° d K°K°, and 10~ 8 for B° d K±K^. CP asymmetry vanishes in the B° d -> K°K° modes, 
because only the penguin operators contribute at leading order. However, if FSI contributes, the above results will 
be changed dramatically. For example, FSI effects could induce large T and P UyC via the rescattering of intermediate 
states DD and irir produced from the tree operators, and Pt via the rescattering of intermediate states KK produced 
from the penguin operators. When IZkk deviates from unity through rescattering processes, the branching ratios 
and CP asymmetries of the B — > KK decays could be enhanced. 

We show how FSI effects modify amplitudes of various topologies in the B — > KK decays in Table I. For more 
allowed intermediate states, refer to fl7|j . It is obvious that the rescattering processes DD(irir) — > KK may be 
important due to the large B — > DD(irir) branching ratios. For example, B(B d — > 7r ± 7r ± ) is of order 10~ 5 . It is then 
possible that FSI effects could be significant enough to increase B(B d — > K^K T ) from order 10 -8 to above 10~ 7 . 
For a similar reason, the rescatering processes could induce large P UyC with the CKM matrix elements V UyC , which, as 
interfered with the penguin contributions, result in sizable CP asymmetry in the B d — > K°K° modes. Hence, large 
CP asymmetry observed in the B d — > K°K° modes and large deviation of the observed B d — » K ± K T branching 
ratios from the PQCD predictions will indicate strong FSI effects. 

III. NONFACTORIZABLE AND ANNIHILATION CONTRIBUTIONS 

PQCD factorization theorem for exclusive nonleptonic B meson decays has been briefly reviewed in [Q. In this 
section we simply sketch the idea of PQCD factorization theorem, concentrating on its application to nonfactorizable 
and annihilation amplitudes in the B — > KK decays. 

In perturbation theory nonperturbative dynamics is reflected by infrared divergences in radiative corrections. These 
infrared divergences can be separated and absorbed into a B meson wave function or a kaon wave function order by 
order [Q] . A formal definition of the meson wave functions as matrix elements of nonlocal operators can be constructed, 
which, if evaluated perturbatively, reproduces the infrared divergences. Certainly, one can not derive a wave function 
in perturbation theory, but parametrizes it as a parton model, which describes how a parton (valence quark, if a 
leading-twist wave function is referred) shares meson momentum. The meson wave functions, characterized by the 
QCD scale Aqcd, must be determined by nonperturbative means, such as lattice gauge theory and QCD sum rules, 
or extracted from experimental data. In the application below small parton transverse momenta fcr are included, and 
the characteristic scale is replaced by 1/6 with b being a variable conjugate to fer- 

After absorbing infrared divergences into the meson wave functions, the remaining part of radiative corrections 
is infrared finite. This part can be evaluated perturbatively as a hard amplitude with six on-shell external quarks, 
four of which correspond to the four-fermion operators and two of which are the spectator quarks of the B or K 
mesons. Note that the b quark carries various momenta, whose distribution is described by the B meson wave 
function introduced above. The six-quark amplitude contains all possible Feynman diagrams, which include both 
factorizable and nonfactorizable tree, penguin and annihilation contributions. A factorizable diagram involves hard 
gluon exchanges among valence quarks of the B meson or of a kaon. A nonfactorizable diagram involves hard gluon 
exchanges bwtcccn the valence quarks of different mesons. That is, the PQCD formalism does not rely on FA. 

The hard amplitude is characterized by the virtuality t of involved internal particles, which is of order Mb, and by 
the W boson mass Myy- The hard scale t reflects the specific dynamics of a decay mode, while Mw serves the scale, 
at which the matching conditions of the effective weak Hamiltonian to the full Hamiltonian are defined. The study 
of the pion form factor has indicated that the choice of t as the maximum of internal particle virtualities minimizes 
next-to-leading-order corrections to hard amplitudes Jl^ j . Large logarithmic corrections are organized by RG methods. 
The results consist of the evolution from Mw down to t described by the Wilson coefficients, the evolution from t to 
1/6, and a Sudakov factor. The two evolutions are governed by different anomalous dimensions, since loop corrections 
associated with spectator quarks contribute, when the energy scale runs to below t. The Sudakov factor suppresses 
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the long-distance contributions from the large b region, and vanishes as b = 1/Aqcd- This suppression guarantees 
the applicability of PQCD to exclusive decays around the energy scale of the B meson mass j|. 

A salient feature of PQCD factorization theorem is the universality of nonperturbative wave functions. Because of 
universality, meson wave functions extracted from some decay modes can be employed to make predictions for other 
modes. We have determined the B and K meson wave functions from the experimental data of the B — > Kit and itit 
decays J?|, and the unitarity angle 4>3 = 90° from the CLEO data of the ratio 19|, 



where B(B^ — > if ± 7r =F ) represents the CP average of the branching ratios B(B^ — > K + n~) and B(B® — > K~n + ). 
It has been emphasized that the B — > hit data can be explained using the same angle cj>3 = 90° in PQCD, contrary 
to the conclusion in |2^j2l[], where an angle larger than 100° must be adopted. In this work we shall predict the 
branching ratios and CP asymmetries of the B — > KK decays in the PQCD formalism employing the above meson 
wave functions and the unitarity angle. 

Factorizable annihilation contributions correspond to the time-like kaon form factor. It is known that annihilation 
contributions from the Oi_4 operators with the [V — A)(V — A) structure vanish because of helicity suppression. 
However, those from the 0$$ operators with the (V — A)(V + A) structure bypass helicity suppression, and turn 
out to be comparible with penguin contributions M. Without FSI in PQCD, strong phases arise from non-pinched 
singularities of quark and gluon propagators in annihilation and nonfactorizable diagrams. Especially, annihilation 



amplitudes are the main source of strong phases [j?]. In the FA and Beneke-Buchalla-Neubert-Sachrajda (BBNS) [£2 23 
approaches, where annihilation diagrams are not taken into account, strong phases come from the Bander-Silverman- 
Soni (BSS) mechanism and from the extraction of the scale dependence from hadronic matrix elements |25|. As 
shown in these sources are in fact next- to- leading-order. As a consequence, CP asymmetries predicted in FA and 
BBNS are smaller than those predicted in PQCD. Nonfactorizable amplitudes have been also considered in the BBNS 
approach, which are, however, treated in a different way. For example, they are real because of some approximation 
in j22j, but complex in PQCD (26). Generally speaking, nonfactorizable contributions are less important compared to 
factorizable ones except in the cases where factorizable contributions are proportional to the small Wilson coefficients 
ai or vanish. 

As stated before, the B^ — » K^K° decays involve both annihilation amplitudes from O5 6 and nonfactorizable 
annihilation amplitudes from 0\ i- Their interference then leads to substantial CP asymmetry in PQCD. The B® — > 
K^RT decays involve only nonfactorizable annihilation amplitudes from tree and penguin operators, such that their 
branching ratios can not be estimated, or are vanishingly small in the FA and BBNS approaches. These quantities 
mark the essential differences among FA, BBNS and PQCD. The comparision of our predictions for the CP asymmetry 
in the B^ — > K^K° decays and for the B^ — > K^K T branching ratios with future data will justify our evaluation of 
annihilation and nonfactorizable contributions, and distinguish the FA, BBNS and PQCD approaches. 

IV. FACTORIZATION FORMULAS 

We present the factorization formulas of the B — > KK decays in this section. The effective Hamiltonian for the 
flavor-changing b — > d transition is given by [E7| 



V2 ^ 

v q—u,c 
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with the CKM matrix elements V q = V* d V q b and the operators 



i=3 



(5) 





= {d l q ) V - 


-A(qjbi)v-A , 


= (diqi) V -A(qjbj) V -A , 


o 3 = 


(dibi) V -A 


J2(q~jqj)v-A , 


O4 = (djbj) V -A y](9j9i)v-A , 
9 


5 = 


(dibi) V -A 


Q 


6 = (djbj)v- A y](qjqi)v +A , 
q 


o 7 = 




A^e q {q 3 qj) V - 
9 


hA j 0$ = -(d i b j ) V - A '^2e q (q : jqi)v+A , 
q 


09- 




A^e q {q 3 qj)v- 


3 _ ^ 

-A i O10 = 7;{dibj)v-A/,eq{qjqi)v-A 



(6) 
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i and j being the color indices. Using the unitarity condition, the CKM matrix elements for the penguin operators 
O3-O10 can also be expressed as V u + V C = ~Vt- The unitarity angle 03 is defined via 



V u b = \V u b\ exp(-£0 3 ) . 
Adopting the Wolfenstein parametrization for the CKM matrix upto 0(A 3 ), 



v ud 






Vcd 


V cs 


v cb 


v td 


V u 


v tb 



1 - ^ 

1 2 

-A 



A A\ 3 (p-ir))' 



we have the parameters 



A: 
A : 
Rb i 



,AX 3 (1 - p-irj) 



0.2196 ± 0.0023 , 
0.819 ±0.035 , 



2 

-AX 2 



AX 2 
1 



(7) 



(8) 



V = 0.41 ±0.07. 



(9) 



For the B^ — > K^K decays, the operators O^ contribute via the annihilation topology, in which the fermion flow 

(c) 

forms two loops as shown in Fig. 1. 0\ 2 do not contribute at leading order of a s . 03-10 contribute via the penguin 
topology with the light quark q — s and via the annihilation topology with q — u, in which the fermion flow forms one 
loop. As evaluating hard amplitudes, an additional minus sign should be associated with the contributions, that 
contain two fermion loops. 01,3,5,7,9 give both factorizable and nonfactorizable (color-suppressed) contributions, while 
02,4,6,8,10 give only factorizable ones because of the color flow. The electroweak penguin contributions from O7-10 
have been included in the same way as those from the QCD penguins 03_6- Obviously, the electroweak penguins are 
less important because of the small electromagnetic coupling. 
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contribute via the annihilation 



The diagrams for the B® — > K ± K T decays are displayed in Fig. 2. The operators 0\ 2 

topology, which contain one fermion loop. 0[ c \ do not contribute at leading order of a s . O3-10 contribute via the 
annihilation topology with the light quark q = s or u, in which the fermion flow forms two loops. In these modes 
02,4,6,8.10 give both factorizable and nonfactorizable contributions, while 01.3,5.7,9 give only factorizable ones because 
of the color flow. Only the operators O3--10 contribute to the — > K°K° modes via the penguin topology with 
the light quark q — s and via the annihilation topology with the light quark q — s or d. The penguin contributions 
contain one fermion loop. The q = s annihilation amplitudes involve two fermion loops, while the q = d annihilation 
amplitudes contain both cases of one fermion loop and of two fermion loops as shown in Fig. 3. 

The B meson momentum in light-cone coordinates is chosen as P\ — (Mg/v2)(l, 1, Ot). Momenta of the two kaons 
are chosen as P 2 = (Mb/\/2)(1, 0, Ot) and P 3 — Pi - P 2 . We shall drop the contributions of order {Mr /Mb) 2 ~ 5%, 
Mjf being the kaon mass. The B meson is at rest with the above parametrization of momenta. We define the 
momenta of light valence quark in the B meson as k\ , where k\ has a plus component , giving the momentum 
fraction x\ — k\ /P^ , and small transverse components ki^- The two light valence quarks in the kaon involved in the 
B — > K transition form factor carry the longitudinal momenta X2P2 and (1 — £2)^2, and small transverse momenta 
k 2 T and — k 2 T, respectively. The two light valence quarks in the other kaon carry the longitudinal momenta X3P3 
and (1 — £3)^3, and small transverse momenta k3*r and — respectively. 

The Sudakov resummations of large logarithmic corrections to the B and K meson wave functions lead to the 
exponentials exp(— Sb), exp(—SK 2 ) and exp(— Sk 3 ), respectively, with the exponents 



S B (t) = s( Xl P+, h) + 2 / %(a s (fl)) 



d\i 



S K2 (t) = s(x 2 P+,b 2 ) + s ((l-x2)P^b 2 ) + 2 I ^ 7 (a.(A)) 

Ji/b 2 A 1 

S K3 (t)=s(x 3 P 3 -,b 3 ) + s{{l-x 3 )P 3 -,b 3 )+2 f ^ 7 Mm)) 

Jl/b 3 A 4 



(10) 



The variable 61, b 2 , and 63 conjugate to the parton transverse momentum kix, k 2 T, and k^T represents the transverse 
extent of the B and K mesons, respectively. The exponent s is written as p9| pl| 



s(Q,b) 



dp 
i/b t 1 



hi I ^) .!«.„(//)) - Binjjm 



(11) 
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where the anomalous dimensions A to two loops and B to one loop are 



7T 



9 3 27 J 3 V 2 



with Cf — 4/3 a color factor, / = 4 the active flavor number, and 7^ the Euler constant. The one-loop expression of 
the running coupling constant, 

47T 

= /W/A 2 QCD ) ' (13) 

is substituted into Eq. ( pd| ) with the coefficient /?o = (33 — 2/)/3. The anomalous dimension 7 = — a s /7r describes 
the RG evolution from t to 1/6. 

The decay rates of B^ — > K^K have the expressions 

r = %^W- (14) 

The decay amplitudes *4 + and A" corresponding to B + — > K + K° and £>~ — > K~K°, respectively, are written as 

^+ = !kv;f^ s) + v;m^ s) + Jbv;f^ u) + v t *M p $ - v:M al , (15) 

A- = f K V t F P{s) + V t M^ s) + f B V t F P(u) + V t M^ ~ VuM a i , (16) 

with the kaon decay constant Jk- The notation F (M) represents factorizablc (nonfactorizable) contributions, where 
the indices a and P(q) denote the annihilation and penguin topologies, respectively, with the q quark pair emitted from 
the electroweak penguins, and the subscripts 1, 4 and 6 label the Wilson coefficients appearing in the factorization 

(u) 

formulas. The nonfactorizable amplitude M. a \ are from the operators 0\ 2 . 

The decay rates of B® — > K^K ^ have the similar expressions with the amplitudes 

A = V t *(M^ + M PW) _ v:Ma2 , (it) 
A = V t (M^ } + M P £) - V u M a2 , (18) 



for £?Jj — > K + K and B® — > K K + , respectively. The notations are similar to those in Eqs. ([15]) and ([Tq). The 
decay amplitudes for B^ — > K a K° and — » K°K° are written as 

A' = f K v; F ^ + v t *M^ + f B v?F%*> + v;(m:$ + M p 3 f + m p J$) , (19) 

A' = f K V t F^ + V t Mlt ] + fBV t F$*> + V t (M P $ + MT + M P J?) , (20) 

respectively. 

The factorizable contributions are written as 

r 46 — r 4 "r r 6 i 

/•I /"OC 

pP{s) _ iQ^CpM^ / dx\dx 3 / 61 (i&i b 3 db 3 (f>B(xi, b\) 
Jo Jo 

x{[(l + X 3 )0K (ar 3 ) + r K {\ - 2x 3 )4>' K {x 3 )\ E$(tW)h e {xi,x 3 ,bi,t> 3 ) 

+2r K ^' K (x 3 )E^(t^)h e (x 3 ,x 1 ,b 3 ,b 1 )} , (21) 

pi pOO 

F P{s) = 32ttC f M| / dxidx 3 bxdb x b 3 db 3 4>B (x u h) 
Jo Jo 

(xi,x 3 ,bi,b 3 ) 

+ [xi<p K (x 3 ) + 2r K {\ - Xl )0' K (x 3 )} E ( e f(tW)h e (x 3 , an, 63, h)} , (22) 

/•l poo 

F P(<,) = 32ttC f M b / dx 2 dx 3 / b 2 db 2 b 3 db 3 
Jo Jo 

xr K {[x 3 cj) K (l - x 2 )<f>' K (l ~ x 3 ) + 2cj>' K (l - x 2 )<j> K {l ~ x 3 )] E { £ ' {t^)h a (x 2 ,x 3 ,b 2 ,b 3 ) 
+ [2<M1 - x 2 )cf>' K (l ~ x 3 ) + x 2 cf>' K (l - x 2 )<I>k(1 - *3)] E { a l ] (t^)h a (x 3 ,x 2 , b 3 , b 2 )} , 

(23) 
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for the light quarks q = u and d. The evolution factors are given by 



Elf(t) = a s (t)a[ q) (t) eMS K2 (t) - S K3 (t)] 



(24) 
(25) 



Notice the arguments 1 — x 2 and 1 — x 3 of the kaon wave functions 4>k and <p' K in Eqs. fl26[ ) and (|23|). The explicit 
expressions of the kaon wave functions will be given in Sec. V, where x represents the momentum fraction of the light 
u or d quark. However, to render the annihilation contributions for q = s and for q = u or d have the same hard 
parts, we have labelled the s quark momentum by x in the latter case, and changed the arguments of the kaon wave 
functions to 1 — x. 

The factorizable annihilation contribution associated with the Wilson coefficient from Fig. 1(c) is identical to 
zero because of helicity suppression as indicated by 



F. 



P(q) _ 



16ttCf M b / dx 2 dx 3 I b 2 db 2 b 3 db 3 
Jo Ja 

x{[-x 3 cf) K (l - x 2 )(f> K (l ~ x 3 ) - 2r%(l + x 3 )4>' K (l - x 2 )(f>' K (l ~ x 3 )] 

xE { a q 4 \t^)h a (x 2 ,x 3 ,b 2 ,b 3 ) 
+ [x 2 (f>K(l ~ x 2 )<p K (l - x 3 ) + 2r%(l + x 2 )</>' K {l - x 2 )<p' K (l - x 3 )] 

xE^\t^)h a (x 3 ,x 2 ,b 3 ,b 2 )} . 



Si) 



(26) 

and the two terms in 



The helicity suppression does not apply to the annihilation contributions associated with a 6 , 
Eq. (23) are constructive. It is easy to confirm these observations by interchaning the integration variables x 2 and 
x 3 in the second terms of Eqs. (|23]) and pq). The factorization formulas for F a i from Fig. 1(a) and for F a2 from 



Fig. 2(a), associated with the Wilson coefficient ai(t a ) and a 2 {t a ), respectively, are the same as F^ q \ i.e., vanish. The 

l 3 



expressions of F^ 3 q ^ from Figs. 2(b), 2(c), 3(c), and 3(d), associated with the Wilson coefficients a 3 q> (t a ) + 
are also the same as F^ q ^ and vanish 



(*«), 



The hard functions h's in Eqs (^l|)-(|2^), are given by 

/l e (xi,X3,6i,fo 3 



K yxix 3 M B bi) 
x [6(h - b 3 )K (^Meh) I Q (^M B b 3 ) 
+6{b 3 - b^Ko {^M B h) I {Vx^Msbi)] 



(27) 



h a (x 2 ,x 3 ,b 2 ,b 3 ) 



(Vx^M B b 2 ) 



6{b 2 
0(b 3 - 



- b 3 )H^> (^M B b 2 ) J {V^M B b 3 ) 
b 2 )H^ (V^M B b 3 ) J {^x~ 3 M B b 2 



(28) 



The derivation of h, from the Fourier transformation of the lowest-order H , is similar to that for the B —* Dir decays 
|p|,P6[. The hard scales t are chosen as the maxima of the virtualities of internal particles involved in b quark decay 
amplitudes, including 1/bf. 

4 1 ) = max(V^Af B , l/b u I/63) , 

4 2) = max{y/xTM B , 1/63) , 
= max{y/x^M B , l/b 2 , 1/63) , 

4 2) = max(^M B , l/b 2 , l/b 3 ) , 



(29) 



which decrease higher-order corrections. The Sudakov factor in Eq. ( |10| ) suppresses long-distance contributions from 
the large b (i.e., large a s (t)) region, and improves the applicability of PQCD to B meson decays. 

For the nonfactorizable amplitudes, the factorization formulas involve the kinematic variables of all the three mesons, 
and the Sudakov exponent is given by S — S B + Sk 2 + Sk 3 - The integration over b 3 can be performed trivially, 
leading to b 3 = b\ or b 3 = b 2 . Their expressions are 
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M P ^=M P ^+M P ^ 



pi pOC 

M P{s) = -Z2kC f ^/2N~ c M 2 b \ [dx] b 1 db 1 b 2 db 2 (f> B (x 1 ,b 1 )4 )K (x 2 ) 

Jo Jo 

x{[{xi - x 2 )4> K {x 3 ) + r K x 3 <f>' K (x 3 )] E { e f(t { < l ) )h^ ) (x 1 ,X2,x 3 , h,b 2 , h) 

+ [(1 -X!-x 2 + x 3 )<j> K (x 3 ) - r K x 3 4>' K {x 3 )] E^\tf ) )hf\x 1 ,x 2 ,x 3 ,b 1 ,b 2 M)} , (30) 

p 1 pOC 

M P{s) = -Z2t:C f ^/2N~ c M 2 b [dx] bid&ifcdMflOsi, J>i)$r(*2) 

Jo Jo 

xr K {[{x\ - x 2 )4> K {x 3 ) + tk{x\ - x 2 - x 3 )<j)' K (x 3 )] 

xE^\t^)h%\x u x 2 ,x 3 MMM) 
+ [(1 -xi- x 2 )(f) K (x 3 ) + r K (l - xi - x 2 + x 3 )(j>' K {x 3 )} 

xE^\tf)hf{ Xl ,x 2 ,x 3 MMM)} , (31) 



M P(q) - M P(q) 4- M P{q) 

JVl ai6 — JVl ai + JVi a6 



M P( A q) = Z2kC f ^2N~ c M 2 b \ [dx] \ hdhhdb^Bixuh) 



f [dx] ( 
o Jo 

x{[x 3 (j) K (l - x 2 )4> K (l - x 3 ) - r\{x\ -x 2 - x 3 )(f>' K (l - x 2 )4>' K (l - x 3 )] 

xE { £'{tf)hf{x u x 2 ,x 3 M,b 2 ,b 2 ) 

- [(xi + x 2 )<f) K (l - x 2 )4> K (l - x 3 ) + r 2 K (2 + xi + x 2 + x 3 )4>' K (l - x 2 )4>' K (l - x 3 )] 

xE { a f(tf ) )hf\x u x 2} x 37 b 1 ,b 2 ,b 2 )} , (32) 

pi pOC 

M P £ q) = Z2-kC f ^2~N c M% \ [dx] hdhhdb^Bixuh) 

Jo Jo 

x{[-r K x 3 (/) K (l - x 2 )<j)' K (l - x 3 ) - r K (x 1 - x 2 )<p' K {l - x 2 )(j) K (l - x 3 )\ 
xEif(tf ) )hf\x 1 ,x 2 ,x 3 ,h,b2,b 2 ) 

- [r K (2 - x 3 )4>k(1 - x 2 )(j)' K {l - x 3 ) - r K (2 - x\ - x 2 )<p' K (l - x 2 )<p K (l - x 3 )] 

xE^'(tf ) )hf ) (x 1 ,x 2 ,x 3 ,b 1 ,b2,b 2 )} , (33) 

with the definition [dx] = dx\dx 2 dx 3 and q — u and d. 
The nonfactorizable amplitudes M P 3 f are written as 



M™=M™+M%* 
<p(q) 



M2 q> = -32nC F ^2N c M' B I [dx] [ b 1 db 1 b 2 db 2 (l) B {x 1 M) 

Jo Jo 

x{[{x\ - x 2 )cj) K (l - x 2 )cj) K (l - x 3 ) + r 2 K (x 1 - x 2 - x 3 )(f>' K (l - x 2 )cf>' K (l - x 3 )] 

+ [x 3 (f> K (l - x 2 )4> K (l - x 3 ) + r 2 K (2 + xi + x 2 + x 3 )(p' K (l - x 2 )4>' K (l - x 3 )] 

xEifitf^hfixtA)}, (34) 

for q = u and d. The evolution factors are given by 

£&>'(*) = a s {t)a { ° ] '{t)eM-S{t)[ h3 = bl ] , (35) 
E { a ?(t) = a s (t)a^'(t) cMS(t)\ b3=b2 ] . (36) 

The expressions of M. a \, -M a2 and M. P 3 q ^ are the same as M. P ^ but with the Wilson coefficients a^itf), a' 2 (tf), and 
a' 3 (tf), respectively. The expressions of M. P 3 ^ are the same as M P 3 f but with the kaon wave functions </>^(l — x^ 
replaced by <j>^ (x^, i = 2 and 3, and with the q quark replaced by the s quark. Notice the difference between the hard 
parts of M P ^ q) and M P ^ q \ which are associated with the O^s operators. The former (latter) corresponds to the 
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fermion flow from the b quark to the d quark in the kaon (the d quark in the B meson), i.e., the case with one fermion 
loop (two fermion loops). That is, the nonfactorizable contributions associated with the structure (V — A)(V + A) 
distinguish these two cases, while those associated with the structure (V — A)(V — A) do not. 
The functions j = 1 and 2, appearing in Eqs. (|3(i|)-(|35|), are written as 

h<f = 



with the variables 



[0{h - b 2 )K (DM B h) I Q (DM B b 2 ) 

+9(b 2 - b x )K Q (DM B b 2 ) h (DMBh) 
xK (D j M B b 2 ) , 



for D) > 



x^H { 1 \J\bJ\M B b 2 ) , for# 2 <0 



ITT 
~2 



0(h - b 2 )H { 1] (FM B h) Jo (FM B b 2 ) 



+9(b 2 - h)H { 1] (FM B b 2 ) J (FM B h) 
xK (FjM B h) , f2 



for F? > , 



H^(J\Fj\M B h) 



for Ff < 



(37) 



(38) 



D = X1X3 , 

D\ = Fl = (xi - x 2 )x 3 , 
D\ = -(1 - x\ - x 2 )x 3 , 
F 2 = x 2 x 3 , 

F 2 = xi + x 2 + (1 — xi — x 2 )x 3 . 
For details of the derivation of hP\ refer to lEql. The hard scales are chosen as 



(39) 



(i) 

d 

(2) 



max(DM B , J\D 2 \M B , l/b u l/b 2 ) , 



max{DM B ,J\D*\M B , l/b u l/b 2 



tf = max(FAf s , J\fJ\M Bi 1/h, l/b 2 ) , 



tf = max(FM B ,\J\F$ \M B , 1/h, l/6a) 
In the above expressions the Wilson coefficients are defined by 

Ci 



ax = C 2 + 
, Ci 

a 2 = Ci + 

, C2 
a 9 = — , 

2 iV c ' 



Nr. 



9i 



M - 



,(?)' _ 



-e a C 



Co 



— Ci 
N c 1 



9^10 



a 



(?) 
4 

(-3)' 



Ci 
1 

N~ c 



9i 

N c 
C3 - 



:e a C< 



q<~"3 1 



Ciq 



9i 

Nr. 



(40) 
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4 9) 


= OH 


C e 
' N c ' 




«$•>' 


1 

= iVc 






4 9) 


= C 8 H 


c 5 




4-)' 


1 


^4 


\e q C 7 



° 7 + K 



(41) 



Both QCD and electroweak penguin contributions have been included as shown in Eq. (41). It is expected that 
electroweak penguin contributions are small, as concluded in [ p2[ . 

The pseudovector and pseudoscalar kaon wave functions <pK and <j>' K are defined by 

'"°Vir(*) = I ^e-^y +1 -(0\u(y+) l5S (0)\«) , (43) 



P 3 ~ — ' J 2VT 



respectively, satisfying the normalization 



1 r l 
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The factor r%, 



<Lr,, k{ .r) : / l ix4>' K { X ) = ^=. (44) 



m 0K M 2 K 

rx = -rr~ > m 0K = ; , (45) 

Mb m s + m,d 

with m s and being the masses of the s and d quarks, respectively, is associated with the normalization of the 
pseudoscalar wave function <j>' K . Note that we have included the intrinsic b dependence for the heavy meson wave 
function (j)B but not for the kaon wave functions 0. As the transverse extent b approaches zero, the B meson 
wave function <pg(x,b) reduces to the standard parton model 4>b{x), i.e., 4>b{%) — 4>B{x,b = 0), which satisfies the 
normalization 

<j> B (x)dx = J f— . (46) 



V. NUMERICAL ANALYSIS 



In the factorization formulas derived in Sec. IV, the Wilson coefficients evolve with the hard scale t that depends 
on the internal kinematic variables X4 and bi. Wilson coefficients at a scale fj, < M\y are related to the corresponding 
ones at \x = M\y through usual RG equations. Since the typical scale t of a hard amplitude is smaller than the b 
quark mass nib = 4.8 GeV, we further evolve the Wilson coefficients from fi = mj, down to /i = t. For the scale t 
below the c quark mass m c = 1.5 GeV, we still employ the evolution function with / = 4, instead of with / = 3, for 
simplicity, since the matching at m c is less essential. Therefore, we set / = 4 in the RG evolution between t and 1/6 
governed by the quark anomalous dimension 7. The explicit expressions of Ci(/i) are referred to Q]. 

For the B meson wave function, we adopt the model (?J 



4> B (x,b) = N B x 2 (l - x) 2 exp 



xM B 



colb 2 



(47) 



with the shape parameter = 0.4 GeV |3^]. The normalization constant Nb — 91.7835 GeV is related to the decay 
constant fg = 190 MeV. The kaon wave functions are chosen as 



10 



4> K (x) = -^=f K x(l - x)[l + 0.51(1 - 2x) + 0.3(5(1 - 2xf - 1)] , (48) 



fa K {x) = -=f KX (l - x) . (49) 



cj>K is derived from QCD sum rules [p4| , where the second term 1 — 2x, rendering <fix a bit asymmetric, corresponds 
to SU(3) symmetry breaking effect. The decay constant fx is set to 160 MeV (in the convention f„ = 130 MeV). 
The wave funcitons 4>b and <f>' K were determined from the data of the B — > Ktt decays . 

We employ G F = 1.16639 x 10~ 5 GeV" 2 , the Wolfenstein parameters A = 0.2196, A = 0.819, and R b = 0.38, the 
unitarity angle fa = 90°, the masses Mb = 5.28 GeV, Mk — 0-49 GeV, and m s = 100 MeV, which correspond to 
moK — 2.22 GeV [Q, and the B d meson lifetime t^o = 1.55 ps (rg- = 1.65 ps) |28[ |. Our predictions for the 

branching ratio of each mode are 

B(B + -> if+if ) = 1.47 x 10 -6 , 
B(B~ -> iV^iV ) = 1.84 x 10~ 6 , 
-> K+R-) = 3.27 x 10~ 8 , 
B(B° d -> K-R+) = 5.90 x 10~ 8 , 
5(5^ -> if°if°) = 1.75 x 10~ 6 , 

B(S2 -> if°if°) = 1.75 x 10~ 6 . (50) 

The above values are lower than those of the B — > tttt decays |B@]. Since the £^ — > K^K^ modes involved only 
nonfactorizable annihilation amplitudes, their branching ratios are much smaller than those of the B^ — » K^K 
and — > K°K° modes. As explained in Sec. II, a large deviation of future experimental data from the predicted 
B d — > K ± K T branching ratios will imply the existence of large FSI effects. 
So far, CLEO gives only the upper bound of the B — > KK decays {HI: 



We also quote the upper bound 



B{B ± -> i^if ) < 5.1 x 10~ 6 , 

-> K±K^) < 2.0 x 10~ 6 . (51) 

B(B° d if iT ) < 1.7 x 10~ 5 , (52) 



from |2^] . Obviously, our preictions are consistent with the above data. 
The CP asymmetries are defined by 



= B(B -> KK) - B(B -> gg) 
CP S(B -> ifif) + B(B -» JfJC) ' 1 ' 

Employing the above set of parameters and fa — 90°, we predict 

A C p(B ± -> i^if ) = 0.11 , 
A CP (B° d -» A-ijfT) = 0.29 , 

A CP (B° -» if°if°) = . (54) 

Basically, the values are of the same order of those in the B — > iT7r decays @. The CP asymmetry in the K°K° modes 
vanishes, because they involve only penguin contributions. Measurements of the CP asymmetry in the B d — > K^K 
can justify the PQCD evaluation of annihilation and nonfactorizable contributions to two-body -B meson decays, and 
distinguish the FA, BBNS and PQCD approaches. The significant CP asymmetry observed in the B d — > K°K° will 
indicate strong FSI effects. 

The dependences of the B — > if if branching ratios on the angle ^3 are displayed in Fig. 4. The branching ratios of 
the K^K° modes increase with fa, while those of the if ± if =F modes decrease with fa. The branching ratios of the 
K K modes are insensitive to the variation of fa. The variation with fa is mainly a consequence of the inteference 
between the penguin contributions and the nonfactorizable annihilation contributions M. a from the tree operators. 
Since M. a \ in Eqs. ( jilf ) and (|l~6| ) and M. a i in Eqs. Ji~7j ) and Jl8| ) contain the Wilson coefficients a[ and a' 2 , respectively, 
which are opposite in sign, the behaviors of the branching ratios with fa in Figs. 4(a) and 4(b) are different. 

The dependences of the CP asymmetries on the angle fa are displayed in Fig. 5. The CP asymmetry in the K°K° 
modes remains vanishing. The CP asymmetry in the K^K^ modes drop suddenly from 70% to zero near the high 
end of fa. Since their branching ratios and the denominator in the definition of Aqp are small, the variation with <f> 
is amplified. Figures 4 and 5 can be employed to determine the range of the angle fa , when compared with future 
data. 
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VI. CONCLUSION 



In this paper we have predicted the branching ratios and the CP asymmetries of all the B — > KK modes using 
PQCD factorization theorem. The unitarity angle 03 = 90° and the universal B and K meson wave functions 
extracted from the data of the B — ► Kir and mr decays have been employed. The dependences of the branching ratios 
and the CP asymmetries on the angle <p3 have been also presented. These predictions can be confronted with future 
experimental data. We believe that these modes can be observed in B factories, which have started their operation 
recently. 

The B — > KK decays are very important for understanding dynamics of nonlcptonic two-body B meson decays, 
such as FSI, annihilation and nonfactorizable effects. In the PQCD formalism FSI effects have been assumed to be 
small. As explained in Sec. II, the B — > KK decays are more sensitive to these effects compared to the B — >■ Kir and 
irir decays. Hence, the comparision of our predictions, especially for the CP asymmetry in the B^ — ► K°K° decays 
and for the B® — > K^K^ branching ratios, with future data provides a justification of the assumption. In PQCD 
the CP asymmetry of the i? ± — > K^K® modes depends on annihilation amplitudes. It has been argued that CP 
asymmetries in the B — > KK decays are small in the FA and BBNS approaches, where annihilation contributions 
have been neglected. Therefore, experimental data of CP asymmetries will distinguish the FA, BBNS and PQCD 
approaches. The B^ — > K modes involve only nonfactorizable annihilation amplitudes, such that their branching 
ratios can not be estimated in FA and BBNS. Future data of these modes can also verify the PQCD evaluation of the 
above contributions. 
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modes 


intermediate 
states 


affected 
topologies 


branching ratios for 
intermediate states j35|,p8| 


data of 

branching ratios 


B + - K+K° 


D+D u 


Pc 


< 6.7 x 10 -3 


< 5.1 x 10~ b 


B» d -> K + K~ 


TY + TY~ 


T,P U 


4.7ii;s ±0.6 x 10_b 


< 2.0 x 10" b 


ttV" 


T,P U 


< 9.3 x 10~° 


K U K U 


Pt 


< 1.7 x 10~ s 


B u d - K U K» 


D + D~ 


Pc 


< 5.9 x 10~ a 


< 1.7 x 10~ b 


7T + 7T~ 


Pu 


4.7+i | ±0.6 x 10~ b 


„u U 

7T 7T 


Pu 


< 9.3 x 10"° 



TABLE I. FSI effects in the B KK decays. 



Figure Captions 



1. Fig. 1: Feynman diagrams for the B ± -> K^K® decays. 

2. Fig. 2: Feynman diagrams for the B d — » K^K T decays. 

3. Fig. 3: Feynman diagrams for the B d — > K°K° decays. 

4. Fig. 4: Dependences of the branching ratios on 03 for (a) the B ± — > K ± K {) modes, (b) the B d — » K^K^ 
modes and (c) the B d — > K°K° modes. The upper (lower) lines correspond to the B (B) meson decays. 

5. Fig. 5: Dependence of CP asymmetries on 3 for (a) the B^ — > K^K modes, (b) the B d — > K^K T modes 
and (c) the B° d -> K°K° modes. 
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Fig.4b 
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